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Abstract. We provide a comprehensive overview of the current theoretical understanding of the dynamic 
structure factor of stiff polymers in semidilute solution based on the wormlike chain (WLC) model. We 
extend previous work by computing exact numerical coefficients and an expression for the dynamic mean 
square displacement (MSD) of a free polymer and compare various common approximations for the hydro- 
dynamic interactions, which need to be treated accurately if one wants to extract quantitative estimates 
for model parameters from experimental data. A recent controversy about the initial slope of the dynamic 
structure factor is thereby resolved. To account for the interactions of the polymer with a surrounding 
(sticky) polymer solution, we analyze an extension of the WLC model, the glassy wormlike chain (GWLC), 
which predicts near power-law and logarithmic long-time tails in the dynamic structure factor. 

PACS. 61. 25. he Polymer solutions - 83.10.Kn Reptation and tube theories - 67.70.pj Polymers 



1 Introduction 

Semiflexible polymers constitute a wide class of biolog- 
ically and technologically important macromolecules. In 
contrast to flexible polymers, their persistence length £p, 
which is the length scale defining the cross-over from a 
locally rodlike to a globally coiled conformation, defines a 
mesoscale much larger than the monomer size. Recently, 
reviving interest in applying scattering techniques to semi- 
flexible polymers results from the need for reliable meth- 
ods for analyzing the mechanical properties of biopoly- 
mers such as actin [1,2], microtubules [3], intermediate fil- 
aments [4] , DNA [5] , fibrin [6, 7] or plant cell wall polysac- 
charides [8], but also peptide fibrils [9], wormlike micelles 
[10], and others. The mechanical properties of biopoly- 
mers in solution have also been studied extensively by a 
number of real space methods such as particle tracking, 
single molecule techniques and microrheology. These pro- 
vide a more intuitive access to the studied molecules than 
the classical scattering methods: even without knowing 
exactly what to expect, one obtains telling pictures and 
videos of the system of interest. Scattering methods, on 
the other hand, have not only the advantage of being ap- 
plicable also to molecules that are too small to be visual- 
ized under the microscope, they moreover allow for a very 
efficient quantitative determination of model parameters, 
with the otherwise laborious ensemble averaging automat- 
ically included. 

Clearly though, their useful application requires a faith- 
ful and accurate representation of the analyzed system 
within a theoretical model that serves to analyze the data. 



The standard model for the theoretical analysis of semi- 
flexible polymers is the wormlike chain (WLC) model, 
which represents the polymer as an inextensible space 
curve with a bending rigidity k = ksTip (in three space 
dimensions). Although this model is generally quite com- 
plicated to solve, essentially all linear and non-linear dy- 
namical properties of fundamental interest - in and out 
of equilibrium - can be obtained via a systematic pertur- 
bation calculation around the weakly bending rod limit 
[11], leaving (in our opinion) little incentive for consider- 
ing models involving uncontrolled approximations. All our 
computations, which will be detailed below, are based on 
the weakly bending limit and all of our results are derived 
in (and require for their quantitative validity) the limit of 
large q£p ^ 1 and qL^ 1. 

While dynamic light scattering (DLS) has proven a 
versatile tool for the study of semiflexible polymer in the 
past [4, 12-15], we feel that the surge of applications and 
the improving experimental accuracy call for a more pre- 
cise evaluation of the theoretical predictions than previ- 
ously accomplished. Though dynamical scaling regimes 
of isolated semiflexible polymers have been known for a 
while [16], a thorough quantitative analysis is presently 
still limited by (1) approximations involved in dealing with 
the hydrodynamic interactions; and (2) serious deviations 
of typical measured structure factors from the idealized 
theoretical predictions for an isolated single polymer. The 
latter are usually due to direct polymer-polymer inter- 
actions, which available models fail to predict quantita- 
tively. The present contribution therefore extends previ- 
ous theoretical work for the dynamic structure factor of 
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stiff polymers in solution by various authors [13, 16-19] 
in several directions. First, we compute via diverse routes 
numerical coefficients that depend on the explicit repre- 
sentation of the solvent hydrodynamics and are needed 
for extracting reliable values for the polymer backbone 
diameter and the persistence length from the scattering 
curves. A new expression for the dynamic mean square 
displacement (MSD) of a free chain fully including the 
effect of hydrodynamic interactions is given. The analy- 
sis sheds some light on the quantitative reliability of the 
method and settles a recent controversy concerning the 
initial slope of the dynamic structure factor [17-19] in fa- 
vor of refs. [17,19]. Further, we extend the analysis of the 
dynamic scattering functions to the practically important 
case that the interaction of the scattering polymer with 
the surrounding polymer solution is not negligible. While 
previous discussions of steric cage and entanglement ef- 
fects were generally based on the tube model [14, 17,20] , we 
introduce a new, superior model, called the glassy worm- 
like chain (GWLC) [21]. Very recently, thanks to experi- 
mental progress, pronounced logarithmic tails due to col- 
lective network dynamics have been discovered in actin 
solutions at low temperature [2] (see also Fig. [2|. They 
cannot be explained within the common tube model and 
have been attributed to a glassy slowdown of the single 
polymer relaxation that, according to the GWLC, results 
from an exponential stretching of the long-time relaxation 
spectrum of an 'ordinary' WLC. The GWLC can poten- 
tially accommodate steric (free volume) as well as sticky 
(enthalpic) [21] interactions between the polymers. As we 
show below, it can account for the observed logarithmic 
tails of the dynamic structure factor. It turns out that the 
presence of the surrounding polymer matrix not only af- 
fects the long-time tails of the structure factor but causes 
severe stretching already at intermediate times. It is there- 
fore indispensable to include the newly derived theoretical 
expressions in any attempt to extract a reliable value for 
the persistence length from scattering data whenever the 
tails of the structure factor exhibit any noticeable devia- 
tions from the ideal stretched exponential form predicted 
for a single free polymer in solution. In practice, this re- 
mark concerns essentially all but the most dilute biopoly- 
mer solutions. 

The remainder of the paper is organized as follows. For 
convenience, the second section contains a brief overview 
of our main results, which are compared to previous pre- 
dictions. The third section summarizes the theoretical model, 
both the WLC and GWLC model, as far as needed for 
the calculation of the dynamic structure factor, detailed 
in section m 



2 Summary of our results 

The coherent dynamic structure factor S{q, t) of a single 
polymer chain is defined as 



Siq,t)/S{q) 



1/e 




Fig. 1. General shape of the dynamic structure factor of a 
single semiflexible polymer in solution. The characteristic time 
T^ = {C±/K)q~'^ separates the initial decay regime with decay 
rate Fg from the stretched exponential tail with decay rate 
Fq , which reflects the structural relaxation. For times t ^ 
Te longer than the entanglement time r^ = {(^±/k,)L'^/-k'^, the 
decay may exhibit a more extreme (logarithmic) stretching due 
to the slowdown of the relaxation of wavelengths longer than 
Le by their interactions with the surrounding (sticky) polymer 
matrix. 



where L is the polymer length, d is the size of a monomer, 
or, more accurately, the backbone diameter, and r(s,i) 
are the monomer positions. The brackets (...) denote an 
ensemble average. Experimentally, the dynamic structure 
factor is obtained from the intensity time-autocorrelation 
function measured in dynamic scattering experiments such 
as DLS with the help of the Siegert relation, or by neutron 
spin-echo spectroscopy. 

One has to distinguish three intermediate regimes of 
the structure factor of a single polymer in solution, as 
schematically indicated in fig. [T] The characteristic time 
scale T^ = {(±/K)q~'^ (Cj. is the transverse friction per 
unit length and q the length of the scattering vector), 
which was already identified in refs. [16-19], separates 
the initial decay regime from the stretched exponential 
regime. While the former is due to the rapid random wig- 
gling of the contour dictated by the thermal forces and 
only weakly (via the hydrodynamic interactions) depen- 
dent on the mechanical properties of the polymer, the lat- 
ter bears the signature of the systematic structural relax- 
ation driven and controlled by the bending rigidity. Also, 
while during the initial decay regime the spatial interfer- 
ence of scattered light from distant monomers along the 
contour contributes, there is no distinction between coher- 
ent and incoherent structure factor at late times [16]. Ul- 
timately, the structure factor depends on time merely im- 
plicitly via the dynamic transverse mean-square displace- 
ment (MSD) 



<5ri(i)^([r^(s,i)-r^(,s,0)]2) 



(2) 



Siq,t) 



1 
Ld 



dsds'(e"i(""(^''*)-""("-'°»), 



of an average bulk monomer (the over-bar denotes a spa- 
^ ' tial average), which we assume to be Gaussian distributed. 
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justified by the equation of motion of a single polymer: 

S{q,t)^S{q,0)exp[-qHri{t)/i] {t:^T^'-). (3) 

Which (and how much) of these different regimes is de- 
tected in a particular experiment depends sensitively on 
q£p. For large q£p ^ 1 the stretched exponential tail dom- 
inates, while it vanishes for q£p « 1 [4]. Finally, under 
practical circumstances, polymer solutions used for scat- 
tering experiments are rarely dilute but rather semidilute 
to achieve a decent scattering intensity Hence, even if the 
scattering wavelength can be made smaller than the mesh 
size, interactions between the polymers usually become 
noticeable for times longer than the entanglement time 
Te, also indicated in fig.fl] 



Short times (i <C t^) For t <^ t^ the dynamic structure 
factor approaches a simple exponential, S{q,t)/S{q,0) ~ 
exp(— r'q t). Below, we present a new discussion of this 
so-called initial decay regime, which yields new results. 
Our derivation of the initial decay rate 



r(o) 

9 



^[c-iogM)] 



(rj « t « T^) (4) 



qualitatively confirms some previous results [17, 19], while 
ruling out others [18]. We also compute numerical values 
for the constant C explicitly via diverse routes and using 
various slightly different forms of hydrodynamic mobility 
tensor, which are summarized in table [l] While the value 
in the lower right corner should be considered the most 
accurate result, the spread of the numerical values serves 
well to demonstrate the sensitivity of the initial decay rate 
to the details of the hydrodynamic modelling, thus sug- 
gesting some reservation against a too naive identification 
of the parameters d extracted from fits to experimental 
data with the physical backbone diameter. 

Taking into account longitudinal contour fluctuations, 
the authors of [18, 19] introduced a new time scale Tq ~ 

{q£p)^'^T^. For times t ^ Tq the initial slope is predicted 
to increase by a factor of two due to contributions from 
longitudinal modes, which are sub-dominant for any fixed 
time t in the limit q£p — > cxd but become important for 
any fixed value of q£p upon taking the limit t —^ 0. Instead 
of repeating the extended discussion, we refer the reader 
to [11, 19] for the longitudinal dynamics. For molecules 
that are reasonably stiff on the length scale q~^ (which 
is a requirement for the application of the weakly bend- 
ing rod approximation on which our whole discussion is 
based) , only the contribution due to transverse motions is 
measurable in practice because of the double scale sepa- 
ration 

/^fVii = {qipy'r^'^^T^ « r^'\^ c. {q£p)-^ « 1 . (5) 

The presence of longitudinal contour fluctuations has to 
be considered in the limit q£p -^ 1, though. 



Long times (i ^ r^) For long delay times the dynamic 
structure factor is given by eq. (l3| [17]. For the WLC, one 
calculates a time-dependent transverse 'correlation length' 
£±{t) = {fit/(^j_y/^ that characterizes the wavelength of 
the modes equilibrating within the time t. The transverse 
friction coefficient per length 



Ci 



Attt] 



-iog[d/£^{t)y 



(6) 



owes its time dependence to the hydrodynamic interac- 
tions that facilitate cooperative relaxations over an effec- 
tive polymer length £±{t) « £±{t). [For a more precise 
deflnition the reader is referred to the discussion following 
eq. ( 30 ) .] The transverse dynamic MSD takes the form 



Srlit)^£lit)/£p^ir^^hr^yq^ 



(7) 



with r^"'' = {kBTq^/^/C±£l^^), resulting in a stretched 
exponential decay of the dynamic structure factor, 

r(l/4) /„,„^ \3/4l 



S{q,t) ex exp 



P (-,'■")■ 



(8) 



in agreement with earlier results [13,17-19]. Equation (Isl 
is a valid expression for times t -^ Tg when the entangle- 
ment constraint is not felt. For simplicity we treated the 
value of a weakly time-dependent parameter arising in the 
approximation to the hydrodynamic interactions as con- 
stant. Whenever the phenomenological parameters need 
to be determined precisely from experimental data, the 
use of the improved expression for the MSD of a WLC of 
finite length given in appendix [A] is compulsory. A com- 
parison of both results is provided in Fig. [31 



Terminal relaxation {t ^ t^) For times longer than the 
'interaction' or entanglement time Tg = {(±/K){Le/Tr)'^, 
defined here as the relaxation time of a mode of (half) 
wavelength equal to the characteristic interaction (or 'en- 
tanglement') length Le, the decay of the structure fac- 
tor slows down even more. This regime, dominated by 
polymer-polymer interactions is observable in semidilute 
solutions with gLg « 1, where the scattering intensity can 
still reasonably be described as an incoherent superposi- 
tion of single-polymer contributions but interactions slow 
down the relaxation at long times. In most applications, 
interactions can in fact hardly be avoided at concentra- 
tions yielding sufficient scattering intensity. As derived in 
detail in section |3.4| and in appendix |B] the GWLC pre- 
dicts a crossover to a near power-law or logarithmic relax- 
ation, in accord with high precision dynamic light scat- 
tering data [2]. This differs from the predictions of the 
dynamic tube model, reproduced in eq. (34 1 below, which 



attributes the slowdown to collisions of the polymer with 
a tube-like viscoelastic cage and predicts an algebraic sat- 
uration of the MSD to a constant plateau value. While not 
unsuccessful in rationalizing experimental data [8,22], the 
tube model lacks the versatility to account for the slanted 
plateaus ubiquitously observed. In contrast, the GWLC 
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Hydrodynamic interactions 


Oseen 


transverse 
mobility 


Oseen + 
diagonal terms 


RP + diagonal 
terms 


equations of motion 


Smoluchowski 


5/6 [17] 


4/3 - 7B 


8/3 - 7B 


- 


Langevin 


4/3 - 7B 


8/3 - 7E 


17/12 + 4/3 - 7i3 



Table 1. Values of the hydrodynamic constant C appearing in the initial decay rate eq. ( 56 1 for different methods of calculation 



and different mobility tensors: original Oseen tensor, Oseen tensor with longitudinal degrees of freedom projected out, Oseen 
tensor with diagonal (Stokes) friction term and Rotne-Prager (RP) mobility function with Stokes friction term. 



attributes the slowing down of the dynamics to an expo- 
nential stretching of the relaxation spectrum characterized 
by a typical free energy barrier height UbTE due to sticky 
interactions with crossing polymers, or ksTT for purely 
steric interactions (cageing/entanglement). We present re- 
sults for the dynamic structure factor for the two limiting 
cases, where only one of these contributions dominates. 
At long times, we write 



5rl{f>T,)^5rl^ 



Xoo) + 5rY'{t) 



(9) 



where 5r\ ^ (oo) — ALl/3£pTT* is the contribution of the 
non-interacting higher bending modes, i.e. the saturated 
static MSD of a WLC of length L^, and the MSD Sr];^{t) 
contains the slow dynamics resulting from the stickiness 
or the steric interactions. With the help of eq. (Is]), the 
dynamic structure factor becomes 



Siq,t)^S{q,0)exp[ 



2r 2 

-1 Srj_,L^_ 



(oo)/4]G£,^(<?,i), (10) 



where 



G£,^(g,t) = exp 



-q 5r 



2.G 



(i)/4 



(11) 



The exponential factor in eq. ([10| plays the role of an 
Edwards-Anderson parameter in the limit of infinitely high 
free energy barriers where G£,j^{q,t) — > 1. 

The GWLC prediction for the MSD of a polymer with 
purely steric interactions is derived in detail in appendix 
[B| and the final asymptotic result is: 



<5rl'«(i) 



AT ^ 1 

^^^{V/'[l + (V-e)exp(^)^] 

- log^/^ [1 + exp{T)T] } {tTe^ /t^ -^ c», J^ < 1). 

(12) 

A result valid for all times is easily obtained with the help 
of eq. ([85]) below. For long times the corresponding tail 
of the dynamic structure factor G';r(g,t) decays as a near 
power-law: 

Gjr{q, t) ex i-"(*) {tTe^ /t^ ^ oo, J^ < 1), (13) 

with a logarithmically time dependent exponent 

q^5r\ ^^{oo) 



a{t) 



4jr3/4 



log-i/^[(i/Te)exp(.F).F]. (14) 



The GWLC prediction for the MSD of a sticky polymer 
at long times is given by eq. (58 1 below. To a good ap- 



proximation valid for short and intermediate times and 



for large £, 



2G 4/13 



^i{~t/TA)+\og{t/TA)] 

(t/TA<f,l<f). (15) 



Here, Ei is the exponential integral function. The most 
salient feature of eq. (15 1 is an intermediate logarithmic 



relaxation that becomes most pronounced for large £ ^ 1, 
which gives rise to a logarithmic decay in the dynamic 
structure factor. 



Ge{q,t)^l 



SlpTT'^ 



JE + log ■ 



t 

TA 



{TA<.t^ £ta) , (16) 



where the arc-length distance A between sticking sites can 
be identified with L^ for strong stickiness. A direct numer- 
ical evaluation and a comparison with experimental data 
show that the logarithmic tail in the structure factor ex- 
tends far beyond the upper bound on t given in eq. ( 16 I 
(compare Figs. [2" 



It turns out that in many cases only one pair of new 
parameters (J-' and L^, or £ and yl) is relevant: for purely 
repulsive interactions A —^ oo, while for (strongly) sticky 
interactions, A —^ L^, and J-' —^ 0. A detailed experimen- 
tal study of the crossover between the two idealized cases 
would be desirable. 



3 The model - equations of motion 

3.1 The wormlike chain (WLC) 

We begin our formal discussion by introducing the WLC 
model and by stating the equations of motion. In the WLC 
model the polymer contour is represented as a continuous 
space curve r(s) with a bending energy 



n 



WLC 






(17) 



The bending rigidity is denoted by k. A key property of 
the WLC is its incxtensibility, expressed by the arc length 
constraint |r'(s)| = 1. This nonlinear constraint renders 
the dynamical equations of motion of the polymer difficult, 
but analytical progress can be made in the weakly bending 
rod limit [11], where the polymer contour is parameterized 
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Fig. 2. Comparison of S{q,t) to the experimentally deter- 
mined dynamic structure factor of (native) F-actin at tempera- 
ture T = 15° C, actin concentration c = 0.7mg/ml. Solid lines: 
DLS data for scattering vectors g[(i.m^^] = 9.62,17.13,29.66 
(from top to bottom) [2] (Original data kindly provided by R. 
Merkel.) Dotted lines: Dynamic structure factor of a GWLC, 
by numerical evaluation of eqs. 
alytical approximations, eqs. ([3l 



By using the an- 
^591 and (IsS), nearly 



indistinguishable fits can be obtained. Values of the fit parame- 
ters: £p = 5.89 M-m (determined for q = 29.66 ^im"^), £ = 43.40 
(g[^m-^] = 9.62), 31.36 (g[^m"^] = 17.13). For the highest q, 
q^A'^/ip > 1. 



by small transverse deflections around the ground state, 
which is a straight line (chosen as the z-axis) . Introducing 
transverse and parallel coordinates. 



r(s,t) = [rj_{s,t),s-rn{s,t)], 



(18) 



the weakly bending approximation is formulated as a per- 
turbation calculation in the small parameter e = r'.., where 
the over-bar denotes a spatial average. The arc length con- 
straint implies rf, w rj^/2. To lowest order, i.e. to order 



.1/2 



there are no longitudinal fluctuations. 



To specify the equations of motion, we need an ex- 
pression for the viscous drag. Hydrodynamic interactions 
are described with the help of the hydrodynamic mobility 
tensor H(r), which relates forces to velocities. We use the 
following Rotne-Prager (RP) form [23]: 



H(r) 



1 



IS{r) 



1 



e{r ^ d){l/r)[l{l + d^ /6r^)+ 
ff{l-d'^/2r^)] (19) 



Here, d is the hydrodynamic diameter. The first term ac- 
counts for the Stokes friction of a monomer. Neglecting 
non-diagonal terms (which only contribute to higher or- 
der in e in the transverse equations of motion), we define 
the following mobility function for the transverse undula- 
tions r±{s,t): 



The linear transverse equation of mot ion including hy- 
drodynamic interactions follows from eq. (17l as [11,17,19] 



dtr±is,t) = f ds' his - s')[-Kr'l'is\t) + ^^(s'.t)], 

(21) 

where ^j_ is the transverse Gaussian white noise. 

The transverse equation of motion is solved by intro- 
ducing normal modes. For simplicity and without serious 
consequences for our results we impose hinged endfjj as 
boundary conditions. 



(22) 



rj_(s,i) = W-^a„(i)sin(fc„s), 



ri=l 



where fc„ ~ m:/L arc the wave numbers. The equation of 
motion eq. (21 1 for / = is rewritten in mode space as 



9ta„(t) = Y^ hnm[-I^Kn^m{t) + lm(^)]: (23) 

nni 

where we have introduced the mobility matrix 

hnm = Y I ^■5^*' sin(fc„s) sin(fc„is')/i(s — s') . 

For n, 771 3> 1 it reduces to hnm = ^nmh{kn), where the 
mode mobility h(k) is the Fourier transform of h(s) [24]. 
For the Rotne-Prager form of the mobility, we thus ob- 
tain the following approximation to the mode friction (for 
fcd< 1): 



a = l/h{k) « 4nTj/[C' - log(fcd)] 



(24) 



The constant C", which takes the value C" = 17/12 — 7^; w 
0.84, is characteristic of the particular form of the hy- 
drodynamic interactions chosen. On this level of descrip- 
tion, the hydrodynamic interactions are completely en- 
coded into this mode dependent friction coefficient of the 
independent normal modes. As usual, the modes relax in- 
dividually and exponentially, 

Cn,n{t) = (a„(t)a,„(0)) = (5„„i(afjexp(-i/T„) . (25) 

Here we defined the relaxation time of mode number n, 
T„ = (Cj^/k)A:~'1. The equilibrium mode amplitudes follow 
from the equipartition theorem as (a^) = 2/^pfc^. 

The most important observable for our discussion of 
the dynamic structure factor is the dynamic part of the 
MSD 

Srl{s,s\t) = 2[{r^{s,0)r^{s\0)) - {r^{s,t)r^{s,0))] 

(26) 



^ For free boundary conditions, we expect some influence of 
the different terminal relaxation times and amplitudes result- 
ing from different values of the low wave numbers fc„ [1] . Our 
discussion is therefore restricted to times shorter than the ter- 



h(s) = S{s)/3TTri + 9{\s\ — d)(l + d /6s )/8TTris. (20) minal relaxation time tl, defined below eq. (271 
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of contour element s with respect to contour element s' . 
In mode space, this takes the form 



6rl{s, s', t) = j-^Yl {cos[kn{s' ~ s)] + cos[fc„(s + s')]} 

P n 

1 - exp(-tn^/rL) 



/c„ 



(27) 



where tl = ti = {C,_i_/ k){L/'k)^. This equation may be 
further simplified by averaging over the variable s with 
s' — s = const. <C L (denoted by an over-bar), which is a 
valid procedure everywhere except in vicinity of size ^j^ (i) 
(to be defined below) of the ends, 



5ri(s,s',t) 



4 



^cos[/c„(s' ~ s)] 



1 



exp(-in''/TL) 

(28) 
In the following we evaluate the above sum on two differ- 
ent levels of approximation. We first present a simplified 
discussion, exactly valid in the limit L ^ oo. In appendix 
\K\ an improved approximation is discussed, and the ex- 
pression given there should be preferred over the follow- 
ing eqs. (29 1, (30) — or their analogue for a WLC of finite 
length, eq. (62) — for quantitative purposes. A compari- 
son of both approximations is shown in Fig. |3] 

We now proceed with the simpler, but less accurate 
expression for an infinite chain. The sum eq. ( p8| is con- 
verted into an integral, and after a change of variables 



w^^-r^ 



i±{t)/d 



\ z{s-s') 
dz cos — — — - — 

1 — exp(— z^) 



(29) 



where we assumed that i± (t) <C L. Here, the weakly vary- 
ing logarithm in the mode friction C-l(^) is treated as a 
(time dependent) constant with respect to k, 



ait) 



Attt] 



C -\og[d/l^{t)] 



(30) 



We have also introduced the two abbreviations l^_{t) and 
Zj^(t) for the (exact and approximate) transverse correla- 
tion length ^j_ (t) = {nt/Ci^f/'^ andZ_L(t) = {nt/'iin^z^Y/'^, 
respectively, with z* being a weakly time-dependent effec- 
tive mode number of order unity. For qualitative purposes 
it is sufficient to substitute J^^{t) -^ q~'^/^(.p in the ar- 



gument of the logarithm of eq. ( 30 1 , which yields a simple 
expression for C-l(*) [19]- 

If the z-integration in eq. (29) is carried out for s = s' , 
the dynamic MSD i5r^(t) = Sr^is.s.t) of a free polymer 
of infinite length is recovered: 



5rl{t) 



£i(t)r(l/4) ^ r(l/4) 
TT^p 3 37rg2 



(r(^)i) 



3/4 



(31) 




Fig. 3. Comparison of different approximations of the MSD 
5r\ yi of a WLC of length A — 1 (im {£p — 1 jxm) with hydro- 
dynamic interactions: numerical evaluation of eq. (70 1 ( — ), 
simple analytical approximation eq. ( 62 1 with time depen- 
dent friction constant (±{t) ( ) and improved approxima- 
tion eq. pi]) (...). 



The full mean square displacement consists of a static 
and the dynamic part, 

Ziri(.s,.s',t)^([r^(s,i)-r^(s',0)]2) 

= (K(s,0)-r^(s',0)]^) + 5ri(s,s',t) 

(32) 

It is easy to see by a Taylor expansion, that in the bulk of 
the polymer the static MSD between two different points 
on the contour is to leading order quadratic in the con- 
tour length, Ar'^{s, s' , 0) ~ e(s — s'Y /2. More precisely, a 
systematic calculation for hinged ends gives [25] 

Arl{s,s')^\l-'L'[{~s'-~sf{l+ 
4[(l/2 - ~sf + (1/2 - ,?)(l/2 - S') + (1/2 - ~s'f] } 

-2|S'-S|3| (S = s/L, S' = s7L). (33) 

The terms in the middle as well as the last one can be 
neglected in the bulk of the polymer, where \s — L/2| ^ 
L/2. For long times, only the dynamic part of the MSD 
evaluated at s = s' contributes to the deacy of the dy- 
namic structure factor, resulting in the incoherent dy- 
namic structure factor [16]. 

In the remainder of this section, we discuss extensions 
of the theory for isolated filaments that address the ef- 
fects of the surrounding solution. For completeness, we 
first briefiy summarize the idea behind and the most im- 
portant prediction of the standard tube model, before we 
describe the alternative GWLC model, which compares 
more favorably to a large body of experimental data. 

3.2 The tube model 

Quasi-static quantities, such as the plateau modulus of 
a solution of semifiexible polymers are satisfactorily ex- 
plained by the tube model, which assumes that the effect 
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of the surrounding network may effectively be represented 
by a tube-like cage, see e.g. [20,26-28]. The tube proper- 
ties are characterized by the entanglement length L^ (the 
arc length over which the bending energy equals the con- 
finement energy) or equivalently by the tube diameter d. 
Both quantities are related by the basic roughness rela- 
tion (P ~ L^el^p '^^ the WLC, which also implies that the 
transverse mean-square deflections from a clamped end 
grow like the third power of the wavelength. 

In the simplest static version of a harmonic tube poten- 
tial, which is added to the Hamiltonian, the relaxation of 
a confined polymer is exponentially suppressed for times 
longer than the entanglement time Tg associated with a 
mode of wavelength k^^ ~ Lg [14]. The saturation of the 
MSD and the dynamic structure factor to their plateau 
values is however too quick if compared to experimental 
observations. This is improved by an extension of the tube 
model that includes dynamic fiuctuations of the tube aris- 
ing from an over-damped homogeneous elastic background 
material [16], which yields 



5rl{t) = (r^) 



l-^Ki)3/^r(-3/4,L^,t) 



(34) 




Here, the prefactors (r^) and (r^^) are the mean squared 
amplitudes of the effective medium and of the polymer, re- 
spectively, uj-i, is approximately (but not strictly) identical 
to the inverse of the entanglement time Te . Theoretically, 
the static MSD, to which the dynamic MSD saturates al- 
gebraically (like i~^/^), is connected to the crossover fre- 
quency uji, via 



(35) 



While this model seems to agree reasonably well with 
experimental data [8, 22, 29] if the amplitudes and the 
crossover frequency are treated as free fit parameters, it 
cannot account for the 'slanted plateaus' generally ob- 
served, which correspond to a very slow terminal relax- 
ation of the MSD. 



3.3 The basic idea of tlie glassy wormlike chain 
(GWLC) 

It is only very recently that high-precision DLS exper- 
iments unambiguously demonstrated a slow logarithmic 
terminal relaxation instead of a plateau in the dynamic 
structure factor of actin solutions at low temperatures [2] . 
This has been interpreted as a signature of a system near 
its glass transition, which only slowly relaxes into equi- 
librium. An interpretation of the logarithmic tails in the 
framework of the established mode coupling theory for 
glasses [30], as suggested by recent simulation studies of 
fiexible polymer blends [31], would require an improba- 
ble fine tuning of parameters into the neighborhood of a 



higher-order mode coupling singularity. It therefore seems 
at variance with the generic nature of the slow relaxation 
in actin solutions, which is found to extend over more than 
three decades in time [2] and observed for a wide range of 
concentrations. For the same reasons, and additionally for 
the lack of any observations of a percolation structure in 
actin solutions, an interpretation in terms of a percolation 
critical point seems unnatural to us. 

In contrast to models attributing the stretching of the 
relaxation spectrum to a critical point, the GWLC [21] 
suggests a very intuitive origin of the slow-down in a rough 
free energy landscape. The basic assumption of the model 
in its simplest form is that short wavelength modes with 
a (half) wavelength A„ = 7r//c„ < A can relax freely, while 
for the relaxation of a mode n with A„ > yl a certain 
number Nn of energy barriers of height EksT have to be 
overcome. This prescription is in the spirit of the generic 
trap models recently favored by many experimental inves- 
tigators of cell mechanics [32,33], but it puts them on a 
more concrete basis. Physically the parameter S is thought 
to arise primarily from direct adhesive polymer interac- 
tions [21]. But in the same vein steric caging and entan- 
glement effects may be cast into the language of free en- 
ergy barriers by a free volume argument. Assigning an Ar- 
rhenius factor for each free energy barrier, the above pre- 
scription gives rise to an anomalous stretching of the WLC 
relaxation spectrum that manifests itself in anomalously 
slow (logarithmic) decay of single-polymer conformational 
correlations at long times. Implementing the correspond- 
ing prescriptions in eq. ( [58| for the dynamic MSD of a 
WLC, one obtains the predictions of the GWLC for the 
late time dynamics. Fig.|4]illustrates the effect of different 
types of interactions of the test chain with its surrounding 
medium on the dynamic MSD, and compare them to the 
predictions of the classical tube model. A similar compari- 
son for the structure factor is presented further below. We 
remark that experimental evidence of logarithmic relax- 
ation (or near-constant loss in the susceptibility spectra) 
observed for synthetic polymer melts [34] could be indica- 
tive of a mechanism akin to that proposed for the GWLC 
in these systems. This observation has recently prompted 
the introduction of the 'glassy Rouse chain' [35]. 



3.4 Direct sticky interactions in the GWLC 

In the simplest version, the GWLC behavior arises from 
very short-ranged (sticky) interactions of crossing poly- 
mers, as they might arise from hydrogen bonds or hy- 
drophobic patches on a biopolymer's backbone, or disper- 
sion forces that are cut off at short distances. If such in- 
teractions are imperfectly screened by electrostatic repul- 
sions, the resulting pair interaction potential will feature 
an energy barrier, which is denoted by EksT [21]. While 
such interactions may under some conditions induce phase 
separation, we concentrate on conditions of moderate at- 
tractions, where they have no thermodynamic but merely 
a kinetic effect via the barrier. In this case, the strat- 
egy of keeping the equilibrium mode amplitudes (a^J of 
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Fig. 4. MSD of a free polymer (— ■ — ), of a polymer in an 
overdamped elastic background (- ••) with (r^) = 0.35 (r^^), 
GWLC with purely steric interactions {■■■), and a GWLC 
with purely sticky interactions according to A = Le, £ = oo 

( ) (equivalent to a static tube) and f = 25 ( — ), all with 

5r\ j^^{oo) — {r\)/3 — ALl/Slpiv'^ = 0.5. The approximate an- 
alytical expressions eqs. ( |15[ ), ( |62[ ) and the asymptotic result 
eqs. (12 1, (62 1 (for t > Tee^'^/JT with A -^ L^) are shown in 
gray. 

the WLC unchanged while modifying its relaxation spec- 
trum should be adequate. (In fact, the question how the 
equilibrium amplitudes are renormalized by the presence 
of a disordered background poses a formidable theoretical 
challenge, but the following analysis suggests that it might 
not be the dominant effect for the dynamics.) 

Whenever such sticky interactions are the dominant ef- 
fect, the relaxation spectrum of the bare WLC is modified 
according to 



r„exp(A''„£') 



n > UA 
n < UA 



L/A 



(36) 



with 



iV„ = (A„/yl - 1) 



(37) 

Here, ua denotes the interaction wave number, so that 
only two new parameters are introduced: the stretching 
parameter or barrier height £ and the interaction length 
A, which has the direct physical interpretation as a typ- 
ical contour distance between sticky contacts. For exam- 
ple, if the polymer contour as a whole is highly sticky, 
A will be identical to the entanglement length L^- When 
the stickiness is only induced at certain places along the 
backbone, e.g. if it is due to an incomplete coverage of the 
backbone with some molecular crosslinker, A may be sub- 
stantially larger than Le ■ The same role is actually played 
by the strength of the attraction in the sticky case, where 
it determines the equilibrium ratio of bound to unbound 
entanglement points and thus A/L^ via a Boltzmann fac- 
tor. In particular in cases where f ^ 1 and/or A ^ Le 
it becomes important to also consider the effect of steric 
interactions between adjacent polymers, known as cage- 
ing or entanglement, which is the subject of the following 
paragraph. 



3.5 Cageing and entanglements in the GWLC 

The GWLC model as introduced above is applicable to 
systems dominated by adhesive interactions. Additional 
contributions to the stretching of the relaxation spectrum 
arise however from steric interactions of the chain with 
its surroundings. Our basic assumption is that for each 
wavelength there are several substantially different confor- 
mations of similar free energy available. In order to relax 
the caged long wavelength modes, the surrounding matrix 
then has to be temporarily pushed out of the way to cre- 
ate enough free volume for the conformational change, a 
process that can again be described by an escape over a 
free energy barrier. More precisely, a wavelength depen- 
dent entanglement volume V\ may be defined as the free 
volume needed by a mode of (half) wavelength A to relax, 
which depends on the magnitude of the transverse excur- 
sions. These may be inferred from the following simple ar- 
gument: With the help of eq. ( 22 1 the MSD of a monomer 



after relaxation of a mode of wavelength A = L/ii is cal- 
culated as 



6r^ 



±,x 



n=L/X 



4L3 



dn 



XVip (38) 



L/A 



(where the spatial averaging has already been carried out). 
Accordingly, we have V\ = XSr'j_ ^ oc A^. The free en- 
ergy barrier, that a mode of (half) wavelength A has to 
overcome to relax in presence of a background polymer 
solution of entanglement length Le is then estimated as 
ksTTX'^/Le- Moreover, it is assumed that modes with 
a transverse MSD smaller than Sr^ ^ can relax freely. 
Similar to the above case with adhesive interactions, this 
leads to an additional stretching of the relaxation spec- 
trum, where the relaxation times of modes of wavelength 
A„ > Le are slowed down according to 



n > He = L/ Le 



T„ exp(iV^J^) n < He 



ith 



N' 



(A^/ie-l)- 



(39) 



(40) 



The dimensionless free energy barrier height J-' is gen- 
erally expected to have a small numerical value, as can 
be seen from the following argument. The interaction free 
energy of a rigid polymer segment with hard core inter- 
actions is estimated as cvkBT/2, where c is the concen- 
tration of collision points and v is the excluded volume of 
the segments. Here we consider semiflexible polymers of 
length Le with purely hard core direct interactions that 
are represented as soft cylinders interacting via an effec- 
tive potential V{r) on a coarse-grained level. Comparing 
their mutual repulsion to a hard core interaction, we esti- 
mate !F ~ 2B2/V, where B2 is the second virial coefficient 
of the solution of soft cylinders. The effective interaction 
potential between the cylinders, is due to a kind of Hel- 
frich repulsion [36,37]. Consider two WLCs of length Le 
that approach each other at orthogonal orientations. If 
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the two axes through the polymer endpoints are held at a 
fixed distance r from each other, their effective interaction 
potential V{r) is found to be [20,38] 



V{r) = -keTlog [erfc(-2%/3r/d)/2 



d = jLl/ip 



From this, the second virial coefficient is calculated as an 
integral over the Mayer /-function. 



B2 



dr 



An 



dn' 

47r 



{exp[-(3V{r,[2,n')]~l}, 

(42) 

where r denotes relative separation between the two axes 
and f2, f2' denote their orientations. The implicit assump- 
tion is made that the potential of the two polymers cross- 
ing at an arbitrary angle is still of the form for orthogo- 
nally crossing polymers. The integral is split into a sur- 
face integral over a parallelogram (with edges of the tube 
length Le) in the plane spanned by two unit vectors in 
the direction of il, fi' (the excluded area) and an integral 
over the distance of the cylinders in a direction trans- 
verse to this plane, and it is assumed that the interaction 
potential only depends on this distance. One arrives at 
B2 = L\d^l%\fi, and using v = 7rLgd/2 for the ex- 
cluded volume of hard cylinders of diameter d and length 
Le [39], we get: 



T: 



1 



/^2V3 



0.2. 



(43) 



For the strength of the steric interactions in Figs. |4][5] 
this value of T is used. It should be noted that T must 
in principle be recalculated if attractive interactions are 
present, in which case its absolute value can be substan- 
tially smaller. For strongly attractive interactions, T may 
turn negative, indicating that the free volume approach 
breaks down. The value for T obtained in eq. (43 1 should 



thus be understood as a rough estimate of an upper limit 
for purely steric interactions, while in the case yl ~ Lg of 
strong stickiness, J^ ^ 0. 



4 The dynamic structure factor 

4.1 Initial decay (i < r^) 

In this section we discuss the initial decay regime of the dy- 
namic structure factor of (isolated) semiflexible filaments. 
Our calculation will be different from that in refs. [17-19], 
which was performed in mode space or Fourier space re- 
spectively, whereas our approach employs a real space rep- 
resentation. 

To perform the thermal average in the dynamic struc- 
ture factor eq. ([T]), the assumption that r_L is a Gaus- 
sian distributed variable is employed, valid for situations 
where the scattering can be traced back to single polymers 
which are described by eq. (21 ), i.e., whenever the scatter- 
ing wavelength is smaller than the mesh size (q^^ ^ ^) 



A 




Fig. 5. Dynamic structure factor of a free polymer (— • — ), 
of a polymer in an overdamped elastic background (- ••) with 
(r^) — 0.35 (r^^), of a GWLC with purely steric interactions 
(■•■), and a GWLC with purely sticky interactions according 

to A — Le and f = 00 ( ) (equivalent to a static tube) and 

£ = 25 (— ), all with 5ri_i^(oo) = {ri)/3 = iLl/Slpn'^ = 0.5. 
The dynamic structure factors obtained from the approximate 
analytical expressions eqs. S, ( |15[ ), (62 1 and the asymptotic 
result eqs. (|3|, ([l2|, ([62| (for t > r^e^^jT with A -^ Le) are 
shown in gray. 



focus on the largest measurable values of q. Neglecting 
longitudinal fluctuations, we get 

S{q,t) ^ ^ J ds J ds'{exp[~qiAri{s,s',t)/A 

+ iq\\{s-s')])o, (44) 
where qj| = qcos{9) and q^ = q^sm'^{9). 



(...)o = (l/2)y d0sin(0) 




(45) 



denotes an average over orientations of the polymer. Tak- 
ing this average, we get in the limit q'^Ar'^{t) ^ 1 (i — > 0) 



S{q,t) ^^Jdsds' j dx[\-q\\-x^)Arl{s,s',t)/A 
-1 

+ ...]exp[igx(s-s')]. (46) 

For i = and to zeroth order in e, the expression for the 
static structure factor is obtained: 

1 
S{q,Q) « ^/d-^ds' Jdxexp[iqx{s-s')] = ^. (47) 



quantitative evaluation of the derived results should thus 



The initial decay rate is defined as 



r(°)^-lhn-log(5(g,i)). 



(48) 



10 



J. Glaser et al.: Dynamic structure factor of a stiff polymer in a glassy solution 



We therefore need to calculate 9fZ\r^(s, s',i). We have 

lim dtArl{s,s',t) = -2lim{rj_{s' ,0)dtrj_{s,t)) 

+ 2lim(r^{s,t)dtr^{s,t)). (49) 



The last term in eq. (49 1 is actually a time derivative 
of an equilibrium correlation function and vanishes con- 
sequently. We replace the first time derivative using the 
equation of motion, eq. (21 1, and thus get an integral over 



an equilibrium correlation function for the transverse co- 
ordinates: 



limdtArl(s,s',t) 

= 2Kliui f ds/i(s-S)(r'["(S,i)r_L(s',0)) (50) 
*^o Jo 

The last correlator is obtained via the equipartition theo- 
rem, which yields 



{vT{S)v^{s'))=2£;'5{S-s'), 



(51) 



takin g in t o ac count two transverse directions. Combining 
eqs. (46), (50 1 and (51 1 we get 



1 
limdtSiq,t) = -^^^ fdsds' fdxh{s-s'){l-x^) 

t — ^0 Zij(X J J 

-1 
X cos(ga;(s — s')). (52) 

With the help of the general identity 

L L L 

Jds J ds' /(s ~s') = j ds f{s){L - \s\) (53) 



-L 



decay rate 



and eq. ( 47 1 we find the general expression for the initial 

(54) 



(0) ^ ^^ d,Siq,t) 
" t-o 5(g,0) 



27ri 



ds I dx{L — \s\)h{s){l — x^)cos{qxs). 

(55) 



Using the mobility function, eq. ( 20 ) , an evaluation of 
eq. (55l for d <C g~^ <C i, an assumption usually fulfilled 



in light scattering experiments on biopolymers, gives 
(0) ksTq^ 



r. 



q.RP 



Qn^rj 



■[C~logiqd)]. 



(56) 



(terms of order (qd)^ and higher have been discarded and 
the limit L — > oo has been taken.) The constant is C = 
4/3-t-17/12-7B = 2.17, where 7^ is Euler's constant. This 
value differs from the previously reported value of C = 5/6 
[17]. We attribute this to our improved treatment of the 



hydrodynamic interactions. However, it also shows that 
the exact value sensitively depends on the details of the 
hydrodynamic model, which has to be taken into account 
for the determination of the hydrodynamic backbone di- 
ameter from measurements of the initial decay rate. This 
is also acknowledged in [19]. While our result confirms the 
previous results found in [17,19], Liverpool and Maggs [18] 
obtained a qualitatively different result, where an addi- 
tional logarithmic factor S{q,t)/S{q,0) — 1 ~ ilog(t) ap- 
pears in the short time behavior of the dynamic structure 
factor, diverging in the limit t —^ 0. This results from their 
too inaccurate treatment of the spherical Bessel functions 
occurring in the discussion in mode space. 

We remark that if one chooses the simpler Oseen ex- 
pression for the mobility function, our result for t he in itial 
decay rate due to transverse displacements, eq. (551, re- 



duces (after carrying out the integration over x and in the 
limit L -^ 00) to a result quoted in refs. [19,24], which 
is derived from the Smoluchowski equation by project- 
ing out longitudinal degrees of freedom from the mobility 
tensor, as proposed also in [17]. This amounts to evalu- 
ating eq. (21) of [19] with their scaling function replaced 
with Ti.{x) = (1 -I- d /dx^) sin(x)/a;. An advantage of our 
approach over the calculation via the Smoluchowski equa- 
tion lies in the fact that the Rotne-Prager corrections can 
easily be implemented in the mobility function. Table [l] 
compares the different predictions for the constant C in 
eq. (|56l). 



4.2 Logarithmic tails {t > Tg) 

Since the dynamic structure factor at long delay times is 
of the (incoherent) form 



S{q,t) ^ S{q,0)exp[^qHrl{t)/4] 



(57) 



its late time behavior is trivially obtained from the dy- 
namic MSD. Upon exponentiating. Fig. His translated into 
Fig.[5J Interestingly, the somewhat limited logarithmic in- 
termediate asymptotics visible in the MSD in the strongly 
sticky limit f ;:^ 1 is thereby extended to much longer 
times, so that a logarithmic intermediate decay seems to 
constitute a quite robust feature of the dynamic struc- 
ture factor of semidilute solutions of sticky semiflexible 
polymers (cf. Fig. pi. The remainder of this paragraph 
is therefore dedicated to a closer examination of the an- 
alytical properties of this intermediate asymptotics. For 
simplicity, we set J- = for the following discussion (the 
case JF > is discussed in appendix BJ). 

We start from the dynamic MSD of a WLC and in- 
troduce the stretched relaxation times f„ eq. ( 36 ) . For 



convenience, we quote the general expression for the MSD 
by converting eq. ( |28[ ) to an integral, after taking the limit 
L,nA ^ 00 with A~^ const.. 



5ri{t) = Srlis,s,t) 






dn 
n 



4 [1 -exp(-VT„)]- 



(58) 
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We perform a change of variables and the dynamic MSD 
is written as 



5ri{t) = 6rlAt) + 5r'r{t) 



where 



4^3 



OO 

|^{l-exp[-(t/r^)n4]}, 



(59) 



(60) 



1 

with TA = {C,^_/ k){A'^ /ti'^) (note that we assume the sim- 
plified treatment of hydrodyn ami c interactions here, for 
an improved discussion of eq. ( 60 1 see appendix \M , and 



5rl^{t) 






dn 



{1 - exp{-(t/r^)n4 exp^f (I/ti - 1)]}}. (61) 



For the slow modes of eq. ( 61 1 , substituting J±^ (t) = A/ir 
in the friction constant Cj-C^is sufficient. The first inte- 
gral can be expressed in terms of an incomplete Gamma 
function, 

4/13 

^rlAi) = 3^ {1 - bit/TAfn-h, {t/TA)]] , (62) 

where b = 3/ 4. ( This term is of the same form as the right 
term in eq. (34 1.) The integral saturates in the limit of 
long times, 

5r\A{^) = 4ylV3^p7r^ (63) 

The second integral is not immediately solvable analyti- 
cally. We approximate it within a certain range of param- 
eters. Writing 

fi = p/4exp[-(£/4)(l/7i - l)]n, (64) 

with the rescaled time t = IJta , we approximatively solve 
the implicit equation for n and get 



(65) 
(66) 



l-(4/f)log(n/7io?i/4) 

4 dn 

^ [1 - (4/£)log(n/noii/4)]2^- 



It is possible to numerically determine the fixed mode 
number tiq. This substitution is valid as long as the loga- 
rithmic term in the denominator of eqs. (651, (66 1 is not 



dominant, i.e. for 5 3> 1. The remaining integral then 
reads: 



6r''f{t) 



4yl3 4 



An , 

^^ [1 — exp(— n ) 
n 



■ log(n/noti/4) 



(67) 



4/13 



{A~B). 



These two integrals can be performed, 

A=- f ^[l- exp(-n4)] ^hjE~ Ei(-i) -I- log(t)] 
t J n t 



(68) 
and 

B=% / 5 [1 - exp(-n4)] log(n/noii/4) 

^ { " (69) 

~ ^[-1/4 - log(7^o)]i + Oitf (i « 1). 

In the limit t/£ ^ 1 we can neglect the contribution of 
B (which for long times is ultimately approximated by a 
logarithmic term in t resulting in corresponding power- 
law-decay in the dynamic structure factor) against A and 
the final result, vali d f or short and intermediate times, 
is provided by eq. (151. Combining eq. (Ill, and 

(t ^ ( X)), w e arrive at eq. 



([151 
observing that Ei(— i) — 
A comparison of the result eq. (151 to the numerically 
evaluated expression eq. (61 1 is shown in Figs. [4| [s] 



As the direct numerical evaluation in Figs. 121 151 shows, 
the logarithmic tail of the structure factor for large £ ex- 
tends well beyond the indicated time regime t <^ £ta- It 
follows from eq. (16 1 that the prefactor of the logarithm 



(the slope of the logarithmic tail in a semilogarithmic plot) 
is inversely proportional to the energy barrier height £, 
which is therefore directly monitored by the slope of the 
tail of the structure factor plotted against Int. Both A and 
£ can be accurately determined from a fit of eq. ([3| to ex- 
perimental data with the 5r\{t) given by the full model 



eqs. (591, (621, (15l 



5 Conclusions and outlook 

We have presented a thorough theoretical discussion of the 
dynamic structure factor of a stiff polymer in a (sticky) 
solution. We claim that our predictions, though in quali- 
tative agreement with previous results, are quantitatively 
superior and provide the missing link to a reliable quan- 
titative analysis of microscopic mechanical parameters of 
polymers (backbone diameter and persistence length) via 
dynamic scattering measurements. Even more important 
seem the prospects of applying quasi-elastic scattering as 
a matchless non-invasive microrheological technique to ex- 
plore with high accuracy the so far poorly understood 
parameter dependencies of the GWLC stretching param- 
eter £. This represents in our opinion one of the most 
promising pathways towards a microscopic modeling of 
the mechanical properties of biopolymer networks and liv- 
ing cells. A particularly rewarding application might be 
the search for an apparently scale dependent persistence 
length of biopolymers [40] , which might arise from the fact 
that these polymers have a fibrous substructure requiring 
a GWLC (rather than a WLC) description already on the 
single polymer level. 
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Appendix A Analytical approximation to 
hydrodynamic interactions 

In this section wc will derive an analytical approxima- 
tion for the mode integral of the MSD of the free GWLC 
modes, eq. ( 60 1 , which is also the part of the MSD which 



Consider the first factor in the integrand. It is 



a3/4(z) = 



-W- 



Az^ 



'0 



3/4 



-7(zo) - log ^ 

,4 \ 1 3/4 



-,3/4 



(z < no) 



(74) 



log 



Az" 



3 7(^0) 



A\og{Az^/nf)_ 
We approximate the first factor in eq. (|74|) for 2 ^ no as 



corresponds to a static tube. We rewrite the integral with 
the help of eq. p4|: 



[-\og{Az^/4)f'^A^IHog^''n, 



SrlAt) 



"0 

4/1'^ f dn 



^_31og(y2z) 



4 log 710 



£p7r4 J n 
1 



/I 
-J {1 - exp[-{t/TA)n:^log{no/n)]} , ^nd the l/log(. . . ) term in the second factor as 



(70) 

where no = ylexp(C")/7rd is an upper mode cut-off cor- 
responding to the finite backbone diameter and ta = 
(47r?7/«;)(yl/7r)'* is the approximate relaxation time of a 
mode of wavelength A. We begin with a substitution of 
variables, z'* = n^log(no/n). The lower bound of the in- 
tegral is then zq = log ' no. With a = 4z''/n^, n{z) is 
determined if a is a solution to the equation: 

^4 



aexp(— a) = 4- 



(71) 



Two solutions for ~a exist: the two real branches of the 
Lambert ly-function. In our case, the dominant contri- 
bution to the integral comes from the mode numbers for 
which n <C no, this corresponds to z > n or a 3> 1. The 
relevant solution is therefore given by a = — IV_i(— 4z'*/no). 
This defines the upper bound of the integral as zi = 
(4e)-i/4no. We then have VF_i(-4zVn^) = -Azyn^, 
and, using the derivative of the Lambert W function, 
WUix) = I^_i(a:)/{.T[1 + W^iix)]}, 



dn = 



1- 



1 



1 



(72) 



The Lambert W function may be asymptotically approx- 
imated by W-i{—x) = 7-1- log(x) for X <C 1 with I7I ^ 
|log(a;)| [41]. Here we choose 7 = W-i{—xo) — log(a;o), 
i.e., we expand W-i{—x) — log(a;) to 0th order around 
xq = —Azq/uq (This means, that for very high mode num- 
bers or very short times the approximation breaks down.) 
For the purpose of numerical evaluation, analytical ap- 
proximations to W-i may be used [42]. Hence we have 
the following expression for the integral eq. ( 70 1 : 



SrlAt) 



AA^ 



-3/4 



dza3/4(^)f 1 



1 



1 — a 



l-exp[-(t/r^) 

^4 



(73) 



-1 



logiAzy^ 



0) 4 log no 
-1 



log(V22) 



41og(no) 



log no 

log(\/2z) 
log no 



(75) 



(76) 



Taken together, one finds, after grouping the terms ac- 
cording to orders of log(z): 

a3/4(2)=43/4log3/4^^ 



[co + ci log z + ©(log z/ log no)^] , (77) 



with 



Co 



^^ 3/ 7(^0) 



4 I 4 log no 



log 2 f^ 3 log 2 



8 log no V 2 log no ^ 

_^ logz ] 
2 log no I 



(78) 



Cl 



4 log no 



1 



7(^0) 



3 log 2 
16 log no V^ log no 



1 



By arguments similar to above, the second factor in the 
integrand of eq. ( 73 1 is approximated for a 3> 1 and z <^ 
no as 



1 



1 



1 



1 



1 — a J A log no 

Consider now the integral 



+ 0(1/ log no)^ 



l-exp[-(VT^)z4] 
dz '-- log z. 



With 



dz 



log z 1-1-3 log Zo 



Q4 



(79) 



(80) 



(81) 
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we can (partially) rewrite the integral eq. (80 1 in terms While the limit exists in terms of a Meijer G function, it 



of an incomplete Gamma function (after a change of vari- can also be performed numerically by inserting very small 



ables): 



1 + 3 log zo 1 
94 16' 



3/4 



dy 



exp(-?/) 



{t/rA)z^ 



y 



7/4 



t 



log{ — ]r[-'S/ A, {t/r a) zt] 



logy 



(82) 



The second term can now be integrated by substituting a 
limit representation for the logarithm: 



exp(-y) f exp(-y)y" 

dy ¥7^— log y = lim / dy 



^7/4 



1 



^7/4 



lim -fr(-3/4+n,a)- 
n-.o n 



(83) 



r(-3/4,a)]. 



values of n, which will be sufficient for our case. With the 



help of eqs. (jeOf, ([TTf, (|79|, (|82), i^, eq. (|73| can now 
be fully approximated and the final result valid for times 
t ^ {ATTri/K)d'^ is given in eq. (84). A comparison of this 
result with the simple approximation given in the main 
text is shown in Fig. [3] 



4/13 



£p7r4 V 4 log no 



Co 
3zl 



1 - 



3z-' 



TA 



r[-3/A, {t/TA)zl] 



-ci- 



l + 31ogzo 1 ("AV jii l{r[_3/4 + n,(t/r^)zo1-r[-3/4,(t/r^)zo1}-logf-)rh3/4,(t/r^)zo1 
16 \taJ l«^o n \ta 



9z3 



(84) 



Appendix B Approximation to free volume 
contributions 



In this section we calculate the MSD of a GWLC with 
purely steric interactions, in the special case J^ > 0, £ — 
0. We may start from eq. ( |58[ ), where we introduce the 
stretched relaxation times t„ for the free volume contri- 
butions according to eq. ( 39 1 . The integral is split into the 



the limit x —f 0, 



dn 1 



j,l/4 



41og^/'^[l + t:rexp(J^-z)] 



1 + t^exp(J^— z) 



dz. (87) 



contributions due to the free high modes Srj^ ^ (t) and due 

2,G/ 



to the slow modes Srj^ (t) , as in eq. ( 59 1 . We have 



Now eq. (85 1 may be rewritten, 



Sr'fit) = 



Sr^-'"-'^ ' 



.7r4 ^3/4 



1 
tpTT J Jl 



(85) 



log* 

dzf{iTe^,z)+ / dzf{iTe^,z) 




= A + B, 



where we substituted A -^ L^- We will approximate this 
integral, by a change of variables, exp(z) = tn'*exp[— J^ 
(1/n^ — 1)], with the rescaled time i = t/Te, which implies 



T/W[iTexp{T-z) 



(86) 



ith 



/(x,z)^ l-cxp[-exp(z)] ^ 



logi/-*[l + a;/exp(z)] 



1 - 



1 



1 + x/exp(z) 



with W being the Lambert function. Using the approxi- 
mation W{x) w log(l -I- x) [42], which becomes exact in 



(89) 

We will in the following restrict the discussion to the case 
J^ < 1 and t 3> exp(—J-')/J-. In this case we may approx- 
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iniate the first term in eq. ( 88 ) by the following integral: 



A: 



^Ll 1 1 

47r4 JF3/4 4 



exp(z) 



log^/''[t>exp(J^-z)] 
1 



tT exp{!F — z) 



(90) 



Here we expanded the double exponential factor to first 
order in exp(z). The integral is done analytically and the 
asymptotic result is 



A 



1 



41og^/'*[t>exp(J^)] 



(tJ^cxp(J^) -^ oo) (91) 



We observe that the double exponential factor in eq. ( 89 1 



may be approximated by a step function, 9{z). The re- 
maining integral of eq. ( 88 ) can then be performed ana- 



lytically, and we arrive at eq. (12 1 of the main text. To 



assess the validity of the approximation made, we first re- 
mark that we require log t 3> 1 . We then compare the two 
terms of eq. (88 1 in the limit <J^exp(jF) -^ cxd and find 



A/B ^ -\og-^[iTcyi^{T)] (i>cxp(J^) ^ c»), (92) 



i.e., the approximation leading to eq. (12 1 is valid for 
log[i>exp(J^)] > 1. 

A comparison of this result to the numerically eval- 
uated expression eq. (851 is shown in Figs. l4| Isj Though 



eq. ( 12 ) is only valid asymptotically, qualitative agreement 



is already reached for t :$> exp(— JF)/J^, T < I. 
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